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Abstract. The mapping between the fermion and spinon compositions of 
eigenstates in the one-dimensional spin- 1/2 XX model on a lattice with TV sites 
is used to describe the spinon interaction from two different perspectives: (i) 
For finite A'^ the energy of all eigenstates is expressed as a function of spinon 
momenta and spinon spins, which, in turn, are solutions of a set of Bethe ansatz 
equations. The latter are the basis of an exact thermodynamic analysis in the 
spinon representation of the XX model, (ii) For N —> oo the energy per site of 
spinon configurations involving any number of spinon orbitals is expressed as a 
function of reduced variables representing momentum, filling, and magnetization 
of each orbital. The spins of spinons in a single orbital are found to be coupled in a 
manner well described by an Ising-like equivalent-neighbor interaction, switching 
from ferromagnetic to antiferromagnetic as the filling exceeds a critical level. 
Comparisons are made with results for the Haldane-Shastry model. 



1. Introduction 

To the few and precious situations, wliere interacting quantum many-body systems 
can be analyzed by rigorous methods in great detail, belong a host of models for 
quantum spin chains. Prominent among them is the XX model for exchange-coupled 
electron spins on a one-dimensional lattice, 



Periodic boundary conditions are assumed throughout. The XX model is the special 
case A = of the XX Z model. The latter also has terms AS'|S'|^]^ in its Hamiltonan. 
The XX Z model is exactly solvable via Bethe ansatz for arbitrary values of A [11 [3- 
The composition of all XX eigenstates by free, spinless Jordan- Wigner fermions 
has been the basis of most advances in the study of the thermodynamics, correlation 
functions, and dynamics [3ll[51El[3[51[ni[ini[IIl[Il[Il[Il[I5]. The ahernative 
composition of the same states provides a framework for the interpretation of dominant 
features observed in the spectra and intensity distributions of dynamical quantities. 
Spinons have spin 1 /2 and semionic exclusion statistics [16j . The one-to-one mapping 
between the fermion and spinon compositions was presented in a recent study of the 
m-spinon dynamic structure factors based on product formulas for transition rates 
[TT] . Here we use the fermion-spinon mapping for different purposes. 
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In preparation of our main themes - interaction and thermodynamics of spinons 
- we introduce alternative quantum numbers for fermion momenta and for spinon 
momenta and spins, along with rules that translate the fermion composition of any 
XX eigenstate into the corresponding spinon composition (Sec. For the energy 
of an arbitrary XX eigenstate, the mapping converts its dependence on the fermion 
quantum numbers into its dependence on the spinon quantum numbers. The resulting 
expression is akin to a coordinate Bethe ansatz (CBA) for the spinon momenta and 
spins (Sec. [3|). From an entirely different perspective, the XX chain is interpreted 
as a set of interacting spinon orbitals with internal degrees of freedom exhibiting 
features akin to electrons in partially filled electronic shells (Sec. [J]). A thermodynamic 
Bethe ansatz (TBA) for spinons is constructed from the CBA discussed in Sec. [3] 
and then analyzed exactly (Sec. [5|). Several aspects of spinon interaction and spinon 
thermodynamics discussed for the XX model are compared with corresponding 
properties of the Haldane-Shastry (HS) model ( [Appendix Ap . 



2. XX spectrum from top dov^rn and bottom up 

The complete spectrum of eigenstates of the XX model is described via 
complementary sets of quasiparticles with different exclusion statistics and with their 
vacua at opposite ends of the energy range: (i) the free Jordan- Wigner lattice fermions 
and the closely related partially interacting magnon solutions of the CBA, (ii) the 
interacting spinons. Only the fermions are free. Therefore, it is desirable to perform 
calculations in the fermion representation. On the other hand, it is desirable to use the 
spinons for the interpretation of the spectrum and the dynamics because the ground 
state (physical vacuum) coincides (for even N) with the spinon vacuum. The magnons 
are important because they are natural products of the CBA, the method on which 
most advances in the analysis of the XX Z model rest [TS] . 



2.1. Magnons and fermions 



The structure of the Bethe ansatz equations (BAE) for the states with magnetization 
M, = N/2 - r of the XX Z model, 



n 
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1, 



undergoes a drastic simplification in the XX limit A — ^ 0: 
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i-iy 



with solutions 



27r 

'n 



(2.1) 



(2.2) 



(2.3) 



provided ki + kj ^ tt holds for all pairs of magnon momenta. The Bethe quantum 
numbers Ii are integers for odd r and half-integers for even r over the range \Ii\ < N/2. 
This simple rule becomes more complicated for A 7^ [TH]. 

In those instances, where magnon pairs with momenta ki + kj ^ n exist, the 
limit A — > in (|2.ip is singular, and nontrivial solutions exist y[9j. The solution of the 
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BAE as a limit process then yields both real and complex-conjugate solutions. The 
critical magnon pairs do not contribute to the energy of the eigenstate, 



£; = ^cosfc,. (2.4) 

i=l 

Thus removing one critical pair from an eigenstate yields a degenerate eigenstate with 
Ak = TT and AMz = 2. This kind of degeneracy is predicted and quantitatively 
described by the sh loop algebra symmetry for the case A = [20l [21] . Critical pairs 
occur only for even N. If N is odd the criticality condition for the associated Bethe 
quantum numbers, li + Ij = N/2 mod N, is impossible to satisfy, given that /i, /2 are 
either both integers or both half-integers. 

The XX Hamiltonian (|l.ip can be mapped onto a system of free, spinless lattice 
fermions by means of the Jordan- Wigner transformation followed by the Fourier 
transform [3]: 

Hxx = ^ cos Pi cl^Cp^, (2.5) 
{p.} 

where the sum is over all sets of distinct fermion momenta from the allowed values 

p-^^fn rn forevenAT, 

N " ' \ {0,2,...,2N -2} for odd Nf ' ^ 

The number of fermions in an eigenstate depends on the magnetization M^: 

M^, < Mz < — 

2 2 ~ 2 

For even N the ground state of TLxx is unique. It has Mz — and thus contains 
Nf = A^/2 fermions. For odd N the ground state is fourfold degenerate. 

How are the fermions and the magnons, which originate from different methods of 
analysis, related to one another? It is evident from (|2.3p and (|2.6p that the non-critical 
magnons have momenta and energies that exactly correspond to fermion momenta 
and energies. Whereas the fermions are free, the non-critical magnons are not. They 
scatter off each other elastically. However, the associated phase shift is tt for all such 
events. This reflects the well-known equivalence between hard-core bosons and free 
fermions in one dimension. 

Regarding eigenstates with critical magnon pairs we note that they are at least 
twofold degenerate within their invariant Hilbert subspace of given quantum numbers 
k and Mz- This degeneracy between states with Afc = and AMz — suggests 
that the magnon eigenbasis for states with critical pairs is related to the fermion 
eigenbasis by an orthogonal transformations within the invariant (fc, M2)-subspaces. 
The critical and non-critical magnons of the XX model can be interpreted as fragments 
of string excitations in the context of the CBA applied to the XX Z model. The limit 
A is singular in the BAE (|2.ip but the singularities can be removed by a basis 
transformation. 

One consequence is that the thermodynamics of magnons as analyzed via TBA 
for the XXZ model becomes equivalent, for A — > 0, to the thermodynamics of free 
fermions Another consequence is that the advances recently reported in the 

calculation of transition rates via algebraic Bethe ansatz for the XXZ model [19] can 
be translated into corresponding advances, for A ^ 0, in the fermion representation 
of the XX model [T7j. 



Nf = -- Mz, --<Mz<-. (2.7) 
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2.2. Spinons 

The unique ground state of Ti-xx for even N coincides with the vacuum for spinon 
quasiparticles. The total number of spinons is confined to the range < Ng < N and 
restricted to be even (odd) if N is even (odd). Knowledge of the numbers of spinons 
with spin up and spin down determines both Ng and Mz. 

N++N^ ^ Ns, N+-N^ ^ 2Mz. (2.8) 

The range of momentum values for spinon orbitals was first determined in the 
context of the HS model [TB] . Minor adaptations are necessary for the XX model [T7] . 
The orbitals available for occupation by spinons are equally spaced at Ak — 27t/N and 
their number is (7V-iVs)/2+l, where A^^ = 0, 2, . . . , iV for even N and Ng ^ 1,3, . . . ,N 
for odd N. The allowed orbital momentum values ar^ 

^, = ^nM, m, = ^,^ + 2,...,N-^. (2.9) 

Each available orbital may be occupied by spinons of either spin orientation without 
further restrictions. The spinon quantum numbers mf are thus integers for even N 
and half-integers for odd N. The XX eigenstate specified by spinon quantum numbers 
{m°'} has wave number 




= ^ mod(27r). (2.10) 



In Ref. [17j we demonstrated how to keep track of the spinons in the fermion 
representation for a system with iV = 4, where the (unique) ground state is the 
spinon vacuum. Here we consider N — 5 for illustration, where the ground state 
consists of four 1-spinon states. The set of allowed fermion momentum states for 
Nf = 0, 1, . . . , 5 is shown in Fig.[TJa) and the sets of allowed spinon momentum states 
for Ng ~ 1,3,5 in Fig. [IJb). An expanded version of Fig. [Ija) is shown in Fig. [5] 
with all 2^ = 32 distinct fermion configurations (circles) grouped according to Nf 
and with the associated spinon configurations (squares) indicated. The V-shaped line 
of Fig. [ija) becomes the forked line in Fig. [2l 

X Note the slight change in convention from Ref. |17| regarding the rrii for odd A^. 




Figure 1. (a) Fermion orbitals available to Nf fermions and (b) spinon orbitals 
available to Ns spinons in XX eigenstates for N = 5. Fermion orbitals can be 
either vacant (open circle) or singly occupied (full circle). The particular fermion 
configuration shown in each row represents one lowest energy states for given 
Nf. Spinon orbitals can be either vacant or occupied by up to Ns spinons with 
arbitrary spin polarization. No specific spinon configuration is shown. 
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The exact spinon configuration is encoded in the fermion configuration as 
described in the foUowing: (i) Consider the V or the fork as dividing the fermion 
momentum space into two domains, the inside and the outside. The outside domain 
wraps around at the extremes {fhi = N mod(A'^) = 0). (ii) Every fermionic hole 
(open circle) inside represents a spin-up spinon (square surrounding open circle) and 
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Figure 2. Fermion configurations of all eigenstates for N = 5. Fermionic particles 
(holes) are denoted by full (open) circles. Spinon particles with spin up (down) 
are denoted by squares around open (full) circles. The fermion quantum numbers 
fhi can be read off the diagram. The spinon quantum numbers m'^ are given 
explicitly and can be inferred from the fermion configuration as explained in the 
text. Also given are the wave number k (in units of 2n/N) and the energy E of 
each XX eigenstate. The spinon momenta kf are discussed in Sec. 13.21 
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every fermionic particle (full circle) outside represents a spin-down spinon (square 
surrounding full circle), (iii) Any number of adjacent spinons in the representation 
of Fig. [5] are in the same orbital of Fig. mb). Two spin- up (spin-down) spinons 
that are separated by £ consecutive fermionic particles (holes) have quantum numbers 
separated by 2£. (iv) The spinon quantum numbers are sorted in increasing order 
from the right-hand prong of the fork toward the left in the inside domain (mf) and 
toward the right with wrap-around through the outside domain (m^). 



3. Interacting spinons 

Naturally, the key to studying the spinon interaction for the XX model is the mapping 
between the fermion composition and the spinon composition of every eigenstate, 
combined with the fact that the fermions are free. Here we study this interaction on 
the level of spinon particles, then, in Sec. IH on the level of spinon orbitals. 



3.1. Energy of spinon configurations 

Consider an arbitrary eigenstate of Tixx for finite (even or odd) N in the fermion 
eigenbasis. Given the configuration of fermion quantum numbers {mi, . . . ,friNf} of 
the state selected, the mapping described in Sec. 12.21 produces a unique configuration 
of spinon quantum numbers with range restricted as in (j2.9p : 

^<m?<m^<...<TO^^ <7V-^, a = ±. (3.1) 

Note that we have sorted them into two groups according to spin orientation and, 
within each group, in ascending order. The number of distinct mf-configurations 
satisfying ((3TT|) for fixed N and is (A^ + l)!/[iVJ(iV + 1 - iVs)!]- Summation over 
Ns yields 2^. All eigenstates are accounted for. Their energies are found to have the 
following dependence on the spinon quantum numbers: 

£;({m+ },{m-_}) = E,{N+,N^) +Y. ^^"^{^I'^D (3-2) 

<T = ± 



with 



vr „ „ TT 

■ TO 



,,„--(A^. + l-2j.), a = ±, (3.3) 
and where the reference energy 

£;o(iV+,A^_)= ^ cos(^toO), (3.4) 

{^°} ^ \^{N + N+-N^ + 2), i(Ar + Af+ - Ar_ + 6), . . . 

...,i(3A^-Af++A^_-2)|, (3.5) 

only depends on A^_^ — i.e. on M^. For even A^, Eq is the energy of the lowest 
eigenstate for given M^, but for odd A^ it does not represent the energy of any 
eigenstate. 

The spinons are not free, notwithstanding the fact that expression p.2p is a sum 
of 1-spinon terms. The spinon interaction is hidden in the sorting criterion p.ip . 
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Moving one spinon into a different orbital will, in general, affect more than just one 
term in p.2p . and switching the spin of one spinon affects all terms including the 
reference energy. 

3. 2. Bethe ansatz equations for spinons 

The very structure of Eq. (|3.2p exhibits features characteristic of a CBA solution 
for spinons undergoing two-body elastic scattering. When we employ the spinon 
quantum numbers (j3.ip in the role of Bethe quantum numbers (BQN) for the spinon 
configuration, Eq. (|3.2p can be rewritten in the form 

E{{m+^},{mJ_}) = Eo{N+,N^) + J2 ^'"^^l (3-6) 

with a universal energy-momentum relation e(fc) = sinfc for spinons, provided the 
spinon momenta fcf satisfy the BAE 

iVAf = ^m^-H ^ -fcf ), z = l,...,7V,, C7 = ±, (3.7) 

o-'=±i=i 

OxxiK - kf) = ^sgn(fcf - kf)d^,,. (3.8) 

Only spinons with parallel spins scatter off each other. All spinon momenta fcf of 
given spin orientation are distinct: 

kl^^{ml-N,-l + 2j,), j, = l,...,7V,, a = ±. (3.9) 

They are sorted in increasing order and bounded as follows: 

^ QiV. - + l) < fc- < fc^^ < ■ • • < /c-^ < ^ (^A^ - iiV, + iV, - . (3.10) 

The distance of any fcf from the upper or lower bound is £(27r/A), where £ = 
0,1,...,{N - Ns)/2 + N„ - 1. The complete sets of mf and kf tor N ^ 5 are 
shown in Fig. [2] 

Even though all spinon momenta of given spin orientation are distinct from one 
another, their exclusion statistics is semionic, not fcrmionic. This is demonstrated by 
applying the defining relation [16j 

Arf, = - ^ WAA^.' (3.11) 

cr' = ± 

for the statistical interaction coefficients to the situation described by (|3.10p . 
taking into account that the number of available momentum states Ad^. with A^^ = 
A'^^ -I- A^_ spinons already present is affected both by the next particle added and 
by the shifting bounds. The result is g^^i =1/2 for all combinations of spinon spin 
orientations. With all these spinon properties in the XX model established we are 
ready to analyze their thermodynamics via TBA from the bottom up. This is the 
theme of Sec. [5] 
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4. Interacting spinon orbitals 

Meanwhile we look at the XX chain from an entirely different perspective. Instead of 
considering individual spinons moving along the chain and scattering off each other as 
described by the BAE (|3.7p we focus on spinon orbitals, specified by orbital momenta 
Ki, and populate them with spinons. 

Each spinon orbital in isolation acquires an energy that varies systematically as 
spinons with spin up or down are added, producing effects not unlike those familiar 
from electronic shells in atomic physics. If two or more spinon orbitals are populated 
with spinons we can express the total energy as a sum of intra-orbital and inter-orbital 
terms. 



4-1- Spinons in one orbital 

There are A^+l states with N spinons. According to (j2.9|) . all N spinons are then in the 
orbital with momentum k = tt/2. These + 1 states all have the same wave number, 
/c = 0, but different values of magnetization, = -N/2, -N/2 + 1, . . . , N/2^ 

The state with all spinon spins up is the fermion vacuum. Flipping one spinon 
spin at a time translates into adding two fermions to the orbitals with the highest 
available energies in the band. The energy levels thus follow from recursive sequences, 
one for even Nf and one for odd Nj: 

(Nf - l)n 

E{Nf) = E{Nf - 2) + 2cos ^ ' , £^(0)^0, £^(1) = 1. (4.1) 

The fermionic particle-hole symmetry implies E{N—Nf) = E{Nf). The level sequence 
resulting from (|4.ip is characterized by spacings ranging between AE ~ 0(1) at the 
bottom {Nf < A^) and AE 0(A^"i) at the top {Nf < N/2). In the limit A^ oo 
we convert Eq. (|4.ip into an integral that yields the reduced energy e ^ E/N as & 
function of the reduced magnetization = M^/N: 

\ , ^ \ 1 , ^ 

e = - cos(7rm^), -- < < +-. (4.2) 



The level distribution implied by (|4.2p can be explained qualitatively by a simple 
microscopic model. The spinon-spin coupling within the given spinon orbital is well 
represented by a ferromagnetic equivalent-neighbor Ising interaction: 

The energy-level spectrum of Jie is 

E, = i Je [A^^ _ (2M,)2] , = y - Nf. (4.4) 

If we set Je = 2/ttN we obtain from (|4.4p the following functional dependence for the 
reduced energy = E^/N on the reduced magnetization m^: 

ee = -(l-4m2), -\<m,<+\. (4.5) 

It shares with (|4.2[) several properties: (i) identical values at = and = ±1/2, 
(ii) a quadratic dependence at \mz\ 1/2, (iii) a linear dependence at — ±1/2. 



§ In the context of the XXZ model at A = 1 all these states are degenerate, forming the multiplet 
with total spin St = N/2. 
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Let us pause here and recall that we started from Tixx j & model of localized spins 
with nearest-neighbor coupling on a ring. The spinon spins, by contrast, are no longer 
localized. Hence their interaction tends to be of long range. For the situation under 
scrutiny here, all spinon spins are almost equivalently coupled, somewhat reminiscent 
of electron spin couplings within atomic shells. 

Now we reduce the number Ns of spinons in the 7r/2-orbital gradually. We 
identify the fermion configurations of this A^s-spinon state for all possible spinon-spin 
combinations. Then we calculate the reduced energy e as a function oi Ug = Ng/N 
and rriz in the limit of iV — > oo. The resulting expression has the form 



1 



cos(7rm2 



2 sin 



1 



(4.6) 



for < Us < 1 and \mz\ < ns/2, in generalization of (|4.2p . The factor in square- 
brackets changes sign at — 1/3. The equivalent-neighbor Ising model for the spinon- 
spin interaction is still applicable. It is merely the effective coupling strength Je that 
now depends on and that switches from ferromagnetic interaction at > 1/3 to 
antiferromagnetic interaction at < 1/3, a Hund's rule of sorts. Inspection of finite- 
N XX chains for which = 1/3 is realized indeed shows that in those particular 
multiplets the spinon-spin coupling is absent. 

For a generic spinon orbital, specified by orbital momentum k, the fermion 
configuration consists of three clusters in momentum space: 



TT 37r TT 
TTWlz < p < ris 



K, 



37r 

T 

37r 
T 



37r 



TT 

2 



K < p < 



StT TT 

K<p< — + -ris 



(4.7) 



Integration of the fermion energy band, cosp, over these three regions yields 

e = — cos(7rm2) 2 sin ^— n^^ sin k — 1 (4-8) 

for < ris < 1, I'TizI < ns/2, and 7rns/2 < k, < tt — 7rns/2. The effective 
equivalent-neighbor Ising exchange constant Je now depends on the filling Us and 
the orbital momentum k. A switch from ferromagnetic spinon spin coupling at low 
filling to antiferromagnetic coupling at higher filling exists only for orbitals with 
7r/4 < K < 37r/4. We have already seen that Je = occurs at — 1/3 (amounting 
to capacity Ug/n'^^^ — 1/3) for the orbital with k = tt/2. As we move away from 
that central orbital to either side, the value of Ug at which Je = occurs decreases 
gradually while the value of Ug/n™'^^ increases. At k = 7r/4 and k = 3tt/4 we have 
Je = at ris = 1/2, which means full capacity (jig/n'^^'^^ — 1). The maximum capacity 
for fixed orbital momentum < k < vr is n™'^^ 1 — |1 — 2k/tt\. 



4-2. Spinons in several orbitals 

The most general spinon configuration involves t orbitals with orbital momenta 
I's < Ki < K2 < • ■ ■ < Kt < TT — i^s and with the spinon content of each orbital 
described by the variables 



ttN, 
2N 



N 



1,2, 



t 



(4.9) 
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The energy expression for this state can be rendered as follows: 

t-i 



^COS^i,i+l |cosi/i^j+i 

i=l 



COS fJ,t.t+l i COSZ^t,t+i 



COS Ki+i — COS K; 



COS Ki — COS Kt 



sm Ki + sm Kt 



sm Ki — sm Ki+i 
1 



}, (4.10) 



where 



j=l j=k j=l j=k 

The first t — 1 terms represent a coupling between nearest-neighbor spinon orbitals in 
orbital momentum space. Each such term depends on the momenta of the two coupled 
orbitals and on the "local" conserved quantities p,i^i+i, I'i^i+i, which play the role of 
coupling constants for nearest-neighbor orbitals. The last term has a slightly different 
structure and depends on the smallest and largest orbital momentum values only. 



5. Thermodynamics of spinons 

In this section we allow for the presence of a magnetic field h in z-direction as 
represented by terms —hSg added to the XX Hamiltonian p.ip or by terms — ^cJ,.Cp. 
added to (|2.5p . The thermodynamic properties of Hxx are derived with least effort 
in the fermion representation [J. From the grand partition function of free fermions, 

Z = [|(l + e-''(™''f*-^)) , l3^{kBT)-\ (5.1) 

we infer the grand potential per site in the limit N oo, 

l3uj{T,n) = - lim iV-ilnZ = -^ rfpln f 1-f e^'^^^'P-''') , (5.2) 

N^oc ZTT J _^ V / 

which translates (for /i = —h) into the Gibbs free energy per site in the spin 
representation, g{T^ h) = uj{T, /i) + fiuf — hniz, where Uf = Nf/N and = 1/2 — n/. 
This result is also obtainable from magnons, namely via TBA applied to the XXZ 
model in the limit A ^ [Tg]. 

Here we demonstrate a different thermodynamic analysis of Ti-xx ■ The results of 
Sec. I3.2l are the basis for an alternative TBA, not from the top down via magnons but 
from the bottom up via spinons. We introduce separate densities in momentum space 
for spinons with spin up and spin down: pa{k) for cr = ±, respectively. For the total 
number of spinons per site we write 

Us^^n^^—^ dkp„{k), (5.3) 

where ^ Nc/N. The integration limits are inferred from (|3.10p : 

fcmm = -C^™^: k^e^x = ^^{1 + (7m^) . (5.4) 

The magnetization (per site) is itself expressible in terms of the Pa{k)'. 



1 V - — 

2^''"'" 47r 



= - > '^"'T = 2~X / dkapa{k). (5.5) 
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Given that all spinon momenta are distinct with allowed values equidistant on a 
prescribed interval, as shown in Sec. 13.21 we express the entropy (per site) in the 
form 



27r ^ 



dk 



p„{k)\np„{k) + (l -p^(fc)) ln(l -p^(fc) 



The internal energy (per site) expressed via the Pa{k) follows from (|3.6|) : 



-y 

27r ^ 



dk Pa{k) sin k cos(7rTO^) 



(5.6) 



(5.7) 



The spinon densities Pcr{k) in thermal equilibrium minimize the Gibbs free energy (per 
site), g — u — Ts — hniz. Solving the variational problem, 5g = 0, for the expression 

/'7r(l + (TT7l2) 



dk 



{pAk) 



sin k + ksT 



Pa (k) In Per (k) 
h 



+ (l - p„{k)J\n(l ~ paikjj - ■^(Tp„(fc)| - i cos(7rm^), (5 



assembled from (|5.3p - (|5.7p . is difficult unless we can remove from the integration 
limits and from the argument in the last term. 

A way out is suggested by the observation that the integration limits of p+{k) 
and p- (k — n) are complements in the Brillouin zone [— tt, tt] . Indeed, if we extend the 
domains (|5.4p of both functions Pa-{k) to the full Brillouin zone via the relation 

p+{k) + p^{k - n) = 1 (5.9) 

we can combine the two integrals in (|5.8p into a single integral either for cr = + or 
cr = — over the entire Brillouin zone with a slightly modified integrand and the last 
term eliminated. Keeping the sum over a we write 

1 



9 



E 



dk 



sin k + knT 



pa{k) In p„{k) 



|pcr(fc) si 

+ (l - p,(fc)) In (l - p^k))] - hapAk)]. (5.10) 

Now the extremum problem is readily solved. The spinon densities in thermal 
equilibrium at temperature T and magnetic field h are 

-1 



Paik) 



_,/3(sin k—ha) 



1 



(5.11) 



Substitution of (|5.1ip into (|5.10p produces the explicit result 

l3g{T,h) = -^ p dfcln(^2cosh(^(sinA:-/i))^ , 

consistent with the result (|5.2p obtained via Jordan- Wigner fermions. Hence the 
thermodynamics of the XX model can be described entirely via spinons. 



(5.12) 
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6. Conclusion and outlook 

It is far from straiglitforward to generalize this study of the spinon interaction and 
spinon thermodynamics to the XX Z model at A 7^ 0. The results presented here 
for the XX limit (A = 0) set the stage for one point from which to attack this 
challenge. A natural second point of attack is the Ising limit (A = 00), where the 
mapping between top-down quasiparticles (ferromagnetic domains) and bottom-up 
quasiparticles (antiferromagnetic domain walls) is again transparent and where the 
latter are again spin-1/2 particles with semionic statistics [22j. 

A promising third point of attack is the Heisenberg limit (A = 1). Its spectral 
properties share key features with those of the HS model owing to common symmetries. 
Given that the XX and HS models, which have very different symmetries, exhibit 
similar degrees of complexity regarding the quasiparticle composition of their spectra 
from top down and from bottom up as well as regarding the spinon interaction 
it is useful to compare several key results established here for the XX model 
with corresponding results known for the HS model. Some relevant HS results are 
summarized in [Appendix A[ Noteworthy similarities and differences are pointed out 
along the way. 



Appendix A. Haldane-Shastry model 

The energy level spectrum of the HS model [231 

Hhs — ^ Jij^i ■ Sj, Jij — J 



N . n{i^j) 
— sm 

TT N 



(A.l) 



on a ring of N sites, is generated from the top down by so-called pseudomomenta that 
represent Yangian multiplet states and satisfy the following set of asymptotic BAE: 

M 

Np, = 27Tl,+TT^sgn{p^-pj), z = l,...,M, (A.2) 

where < M < [7V/2] (integer part of N/2). The BQN li are integers for odd M and 
half-integers for even M on the interval 



^(M + l)</, <iV-i(M + l). 
The solutions of (jA.2[) are of the form 

Pi^^rhi, mi e {1,2, . . . , - 1}, m^+i 
The wave numbers and energies of HS levels are 

M 



ruj > 2. 



(A.3) 
(A.4) 



27r 



M 



^^^Y."^^ mod(2^), S = 2(^j^m,(m, -A^), (A.5) 

i=l i=l 

where Vs = tt J/2 and where the origin of the energy scale is set to coincide with the 
vacuum of pseudomomenta (M = 0). The lowest energy level contains the maximum 
number of pseudomomenta. This level is unique if N is even and fourfold degenerate 
if N is odd. The pseudomomenta play a role similar to the fermions in the XX model. 
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(A.9) 



From the bottom up the HS spectrum is generated by spinons. The foUowing 
relations hold between the number of spinons, A^^, the number of pseudomomenta, 
M, and the magnetization, Af^: 

N++N^^Ns = N-2M, N+-N^ = 2M^. (A.6) 

The number of spinons is determined by the number of pseudomomenta alone. In 
Ti-xx the number of fermions alone does not determine the number of spinons. The 
unique HS ground state for even A^ (spinon vacuum) has energy 

The number of orbitals available for occupation by spinons is (A — As)/2 + 1, where 
A^s is even (odd) for even (odd) A^. The allowed spinon orbital momenta are 

t< = -^m,, m^ = ~M,~~M + 2,...,+M. (A.8) 

A generic HS eigenstate may be specified as follows: 

-i(A^ - A,) < mi < ma < • • • < tun, < ^{N - A^,) 

(Tl, (72, . . . , (TNs, O-j = ± 

The use of separate sets of quantum numbers for spinon momenta (m^) and spinon 
spins (ct.;) is the natural choice for the HS model. In the XX model we used separate 
sets of momentum quantum numbers for spin-up spinons and spin-down spinons. The 
wave number and the energy of the HS eigenstate (|A.9p are independent of the (7^ and 
depend on the rm as follows [121 US] : 

fc = Af7r+-^^mj mod(27r), (A.IO) 

i=l 

1 

E-E,, = Em + e(m,) + ^ ^ V{m, - m^). (A.ll) 

2=1 i<j 

The first term in (jA.lip depends (for given A^) only on the number of spinons present, 

TTV ( 1 1 1 

Em = ^ I -^N{N^ + 2) - -M [3N{N - 1) - Ul^ -h 6M + 4] L (A.12) 

the second term depends on the orbital momenta, 

e(m.) = ^ (M^ - m,^) , (A.13) 
and the third term describes a pair interaction of sorts, 

— V (raj — m, = — ^ 

A representation of Yangian multiplets that describes the pseudomomcntum 
content and the spinon content simultaneously is the motif as illustrated in Table ED 
for N = 5. The motif of an HS eigenstate consists of binary strings of length A^. 
The elements of each permissible string are a '10' (pseudomomcntum) and a '0' 
(spinon). All consecutive 'O's that do not belong to a '10' represent spinons in the same 
momentum state. Consecutive 'lO's represent pseudomomenta with Am^ = 2. Every 
'0' between two 'lO's increases Am^ by one unit. Pseudomomenta with increasing 
fhi are encoded by successive 'lO's read from left to right. Spinons in orbitals with 



;V{m, - m,) = -l (M - |m, - m,\) . (A.14) 
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Table Al. Motif, pseudomomentum quantum numbers fhi, spinon orbital 
momentum quantum numbers mi, wave number k (in units of 2n/N), energy 
E — Egy (in units of TrVg/N^), spin content, degeneracy, and spinon momenta ki 
(in units of n/N) of all Yangian multiplets for N = 5. 



motif 


mi 


mi 


k 


E — Esv 


spin 


deg. 


^1 


010 10 


2,4 


2 


1 


5 
2 




1 

2 




2 


2 


10 010 


1,4 








13 
2 




1 
2 




2 





1010 


1,3 


-2 


4 


5 
2 




1 

2 




2 


-2 


0010 


3 


-1,1,1 


3 


21 
2 




31 = |€ 




6 


-3,1,3 


010 00 


2 


-1,-1,1 


2 


21 
2 




3i = |e 




6 


-3,-1,3 


000 10 


4 


1,1,1 


4 


29 
2 




3 
2 




4 


-1,1,3 


10 000 


1 


-1,-1,-1 


1 


29 
2 




3 
2 




4 


-3,-1,1 


00000 




0,0,0,0,0 





45 
2 




5 
2 




6 


-4,-2,0,2,4 



increasing are encoded by successive 'O's (separated by at least one '10') read from 
right to left. 

The spin content of any given Yangian multiplct can be read off the binary motif 
by recognizing the multiplets of the quantum number St representing the total spin 
in each spinon orbital. In the XX case the motif pertains to individual eigenstates 
(Fig. [2]) and encodes a specific spinon spin configuration. 

Just as in the XX model, the energy expression (jA.lip rewritten in the form 

E-E,,^ Eo{Ns) -^Yl "^»("^» -Ns-l + 2i), (A.15) 

1=1 

EoiNs) = (^) ^N, (3iV^ + 4 - n!) . (A.16) 

is suggestive of a CBA for spinons. The with range (jA.9[) become the BQN and 
(jA.lSp turns into 



E-E,, = Eo{N+ + iV_) - ^ V Kjk, , (A.17) 

TT ^ — ' 



i=i 



where the spinon momenta 



k, = ^{m,- Ns'l + 2t) (A. 18) 
are the solutions of the BAE: 

Nh=TTmi+Y^0Hs{ki-kj), i = l,...,Ns, (A.19) 

Onsih - fcj) = 7rsgn(fc, ~ k^). (A.20) 

Two differences from the XX model, Eqs. (I3.7|) - (|3.8|) . are noteworthy: (i) we are 
dealing with just one set of BQN; (ii) all spinons scatter off each other, not just those 
with the same spin orientation. 
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The HS model too can be interpreted as a set of spinon orbitals, each specified 
by an orbital momentum . Each such orbital is filled with spinons of arbitrary spin 
orientation up to a certain capacity. The available orbitals again depend on the total 
number of spinons present. The energy expression for any Yangian multiplet only 
depends on the orbital momenta and fillings. Conclusions can again be drawn about 
the energetics and interaction of spinon orbitals. 

Beginning with the case where all spinons are in the same orbital, we express the 
reduced energy, e = (i? — Esv)/ (ttvsN), as a function of rig = Ng/N and k = K/27r in 
the limit N —^ oo: 

e - ^n, (3 - - 48^^) , (A.21) 

for comparison with the XX result (|4.8p . which also depends on the magnetization. 
The generalization of the HS expression (|A.21[) to t spinon orbitals is again structurally 
similar to the XX result (|4.10p . except for the absent spinon-spin dependence: 

e = ^ - [{ai,i+i - Rif - (aj,i+i - Ri+i)^] + 2 [{a^^i+l - k^+i)^ - {ui^i+i - k^)^] 



3 

=1 



- [{at,t+i - Kt)^ - (ao,i - Kif] + 2 [(ao,i - Kif - {af,t+i - Kff] , (A.22) 



where = M^'^/{2N), 



i t 

az,k^--7:^Vj + i:^V]^ i = 0,l,2,...,t, (A.23) 



4 2^ 2^ 

j=l j=k 



and ivi"'^ is the number of spinons (with arbitrary spin orientation) in the orbital with 
reduced momentum Rj. 

The thermodynamics of spinons for Tins was reported by Haldane along two 
different paths. The approach taken in Ref. [26j uses the spinon orbital momenta 
Ki (named fc) as the independent variables. Since there are no restrictions on the 
occupation of available spinon orbitals for given Ns , a bosonic version of the entropy 
functional is used. The approach analogous to the one taken in Sec. [5] for Hxx would 
use the spinon momenta ki as the independent variables. Since all ki are distinct, a 
fermionic version of the the entropy functional would have to be used as in (|5.6p . The 
approach taken in Ref. 27J introduces rapidities, Xi = {2mi + 2i — Ns)/N, which again 
are all distinct, but have the advantage of being confined to an interval with limits 
that are independent of Ng. Haldane's result for the Gibbs free energy per site. 




sinh(f [1 + M^)]^ 
sinh ^ 



(A.24) 



sinh (^^K^)^ = sinh ^e-''^°(^), eo(x) = ^7:v,{l ~ x^). (A.25) 

is to be compared with the corresponding XX result (|5.12p . The distributions of 
spinon rapidities inferred from (jA.24[) . 



tanh 



^-^^^^ 2 1-e-W+.4 ' I'l^l' (^-26) 
are to be compared with the distribution of spinon momenta (j5.1ip in the XX model. 
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